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Abstract 

We consider the model for the distribution of a long homopolymer in a potential 
field. The typical shape of the polymer depends on the temperature parameter. 
We show that at a critical value of the temperature the transition occurs from a 
globular to an extended phase. For various values of the temperature, including 
those at or near the critical value, we consider the limiting behavior of the polymer 
when its size tends to infinity. 
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1 Introduction 

The goal of this paper is to analyze various critical phenomena for a model of long ho- 
mogeneous polymer chains in an attracting potential field. The model exhibited here 
demonstrates a phase transition from a densely packed globular phase at low temper- 
atures to an extended phase at higher temperatures. In the latter phase, the thermal 
fluctuations overcome the attraction between monomers and the chain takes on the shape 
of a 3(i random walk or Brownian motion with a typical scale 0{-\/T) where T is the 
length of the polymer. A real life example of this phenomenon is that of albumen (egg 
white). We describe a rough picture of this situation. The physical reality is more com- 
plex as there are present several types of protein with different critical points. However 
in a simplified version, at room temperature the albumen is in the globular state and as 
a result, it forms a viscous, translucent liquid. However, at higher temperatures (around 
60 — 65° C) there is a transition of the albumen to a diffusive (extended) state resulting in 
an opaque semi-solid material. While this transition may be reversible for an individual 
polymer, in the aggregate, the polymer strands in the diffusive state become interwoven 
and form chemical bonds with each other and can not return to the globular state when 
the temperature is decreased. 
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It is worthwhile recalhng Gibbs' philosophy of phase transitions. Start with a system 
of finite size T. The configuration space S-p = {x(-)} denotes all possible states x(-) of 
the system. The space St is equipped with a reference measure Po,t which corresponds 
to infinite absolute temperature (in our case, the inverse temperature (3 = 0). The 
configurations satisfy boundary conditions which reflect the interaction of the finite system 
with its environment. This system is endowed with a Hamiltonian Ht giving the energy 
Ht{x) of the state x. For /3 > 0, the Gibbs measure P/3^t is given by the density 

^P^^ ^ ^ _ exp{-l3HT{x)) 

where 

Z0,T= I exp{-(3HTix))dPo,T. (2) 

When T < CX3, the measure P/3,t and the thermodynamic quantities associated to P^j^t 
are analytic functions of (3. 

Now let T — i> oo. In typical situations, there is a critical value Per such that for /? > Per, 
there exists a unique limiting measure on E, the space of infinite configurations, and 
this limiting measure is independent of the boundary conditions on Ey. Moreover, P^ and 
its relevant thermodynamic quantities are still analytic functions of P for P > Per- One 
manifestation of the phase transition is the non-uniqueness for P < Per of the limiting 
measure as T ^ cxd as it has dependence on the boundary conditions on Sy. Another is 
the non-analyticity of thermodynamic quantities associated to P/j as a function of p. The 
mathematical characterization of the phase transition in terms of non-uniqueness of the 
limiting Gibbs measure traces its history to the works of Dobrushin p| and Ruelle [7] . 

Modern physical theories predict that near the critical point P = Per the limiting 
Gibbs measure P/j must be invariant with respect to renormalizations of the system (self- 
similarity). This idea is related to the two-parametric scaling by Fisher [3j for P near 
Per- Another important fact is that critical behavior as P ^ Per of the physical system 
demonstrates universality, that is the same behavior holds for a wide class of Hamiltonians. 

The most essential part of the present paper is the detailed description of the polymer 
chain near the critical point and the establishment of the physical ideas of universality 
and self-similarity for our particular model of homopolymers. 

2 Description of the Model and Results 

A continuous function x : [0, T] M^, x{0) = 0, will be thought of as a realization of the 
polymer. The parameter t G [0,T] can be intuitively understood as the length along the 
polymer (although the functions x = x{t) are not differentiable and the genuine notion of 
length can not be defined). 
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We assume that for /3 = 0, the polymer is distributed according to the Wiener measure 
Pqt on T,T = C([0, T], M''). For an infinitely smooth compactly supported potential 
V G C^(R'^) and a coupling constant /3 > 0, the polymer is distributed according to the 
Gibbs measure P/3,t, whose density with respect to Po,t is 

dtr'o,T Zf3,T 

In other words, the Hamiltonian Ht is given by Ht = — Jq v{x{s))ds. The normalizing 
factor ^/3,T, called the partition function, is given by 

Zf3,T = I exp(/? / v{x{t))dt)dVo,T{^) = Eo,Te-^^^. (3) 

Jc([0,T],K'') Jo 

It will be usually assumed that the potential is nonnegative and not identically equal to 
zero. We shall be interested in the prevalent behavior of the polymer with respect to the 
measure P/j^t as T — oo. 

We shall see that there are two qualitatively different cases corresponding to different 
values of (5. Namely, for all sufficiently large values of f3 there is a limiting distribution 
for x{T) with respect to P/3,t- Moreover, for each positive constant s and each function 
S{T) such that S{T) oo and T — S{T) ^ oo as T — > oo, the family of processes 
x{S{T) + t), t E [0, s], with respect to either measure P/3,t or P/3^r(-|x(r) = 0), converges 
to a Markov process as T — oo. The generator of the limiting Markov process and its 
invariant measure are written out explicitly in Theorem 18.31 Since x{S{T)) and x{T) 
converge to limiting distributions and thus typically remain bounded as T — ^ oo, we shall 
say that the polymer is in the globular state. 

If /3 > is sufficiently small and d > 3, then the family of processes x{tT)/\/T, 
< t < 1, defined on (C([0, T], M"^), P^ t), converges to a Brownian motion on the 
interval [0,1] (Theorem 19.21) . In this case we shall say that the polymer is in the dif- 
fusive state. Similarly, the family of processes x{tT)/\/T, < t < 1, defined on 
{{C{[0,T],M.'^),Pj3^T{-\x{T) = 0)), converges to a Brownian bridge on the interval [0, 1]. 

We shall see that there is a number f3cr (called the critical value of the coupling 
constant) such that the polymer is in the diffusive state for f3 < j3cr and in the globular 
state for P > P^r- The value of Per and the behavior of the polymer when P is near P^r 
depend on the dimension d and on the potential. In particular, we shall see that Per = 
for d = 1,2 and Per > for (i > 3. 

Of particular interest is the behavior of the polymer when P = Per- In this case the 
appropriate scaling is the same as in the diffusive case, that is we study the family of 
processes x{tT)/\/T, < t < 1. We shall find the limit of this family as T — > oo. It 
turns out to be a Markov process with a non-Gaussian, spherically symmetric transition 
function (Theorem 110.61) . The transition function of the limiting Markov process will be 
written out explicitly. 
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In order to determine whether the polymer is in the globular or diffusive state for a 
given /3, we shall look at the rate of growth of the partition function ^/3,t- Namely, let 



MP) = lim ^S^. 

It will be demonstrated that the limit exists and is equal to the supremum of the spectrum 
of the operator Hp = + (3v : L'^{W^) L^(R°'). The infimum of the set of (3 for which 
^o{P) > is equal to Per- It will be seen that Xo{Pcr) = is an eigenvalue of Hp^^ in 
dimensions d> 5, and corresponds to a ground state of Hp^^ in dimensions d = 3,4. 
The paper is organized as follows. 

In Section [3] we consider finite T and show that {x{t),0 < t < T} is a time- 
inhomogeneous Markov process with respect to the measures P/s^t and P^,t {-IxiT) = 0). 

In Section m we prove the existence of the critical value of the coupling constant. In 
Section[5]we analyze the properties of the resolvent of the operator Hjj which, in particular, 
will be needed to study the asymptotic properties of the partition function. 

In Section Owe shall examine the asymptotics of Ao(/?) when j3 I Per and show it has 
the following asymptotic behavior as /3 | Per, 

CiP - Perf, d = 3, 
XoiP) ~ { C,{P - Per)/ HI/ {P - Per)), d = 4, 
CdiP-Pcr), d>5. 

These asymptotics demonstrate universality in that they depend only on dimension. The 
constants Cd, d > 3, are not universal however. In Section [7] we find the asymptotics, as 
T — i> oo, of Zfs^T- In particular, when P > Per, we shall find that 

j^^q\o{P)t ^qj, some 

constant kp, while for P < Per, Z^^t has a finite limit as T oo. Finally, when P = Per, it 
turns out that Z^t ~ k^T^^"^ for d = 3, Z^^t ~ k4^T/\nT for d = 4, while Z^^t ~ k^T for 
d> 5. We also give asymptotics of the solutions to the parabolic equation du/dt = H^u. 

In Sections [H [9]and[T0l we describe the behavior of the polymer for P > Per, P < Per 
and P = Per, respectively, establishing the convergence results mentioned above. 

Some of the results presented above have been obtained by Cranston and Molchanov 
in pp for the discrete model with the potential concentrated at one point. The analysis 
was based on explicit formulas for the solution of the parabolic equation with such a 
potential. The current results demonstrate that the behavior of the polymer is "universal" 
with respect to the choice of the potential. Another essential feature of this paper is the 
detailed analysis of the behavior of the polymer when P = Per- We refer the reader to the 
review of Lifschitz, Grosberg and Khokhlov |5] for a wealth of information and ideas on 
polymer chains. 
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3 Time-inhomogeneous Markov Property 

First we define as the fundamental solution of the heat equation 

^(t, y, x) =^A^p0{t, y, x) + (5v{x)pp{t, y, x), 
Pf3{Q,y,x) =S{x-y). 

In this section we shall prove that with respect to the measure -P/3,t, the process {x{t), < 
t < T} is a time-inhomogeneous Markov process. Since we shall point out the link be- 
tween non-uniqueness of Gibbs measures and phase transitions it will be necessary to also 
consider the transition mechanism for the process {x{t),0 < t < T} under the conditional 
measure Pp^T (■^(T) = 0). Namely, we will show that the free boundary condition corre- 
sponding to the measure P^ t and the pinned boundary condition corresponding to the 
measure P^ ^- (■|a;(T) = 0) lead to different Gibbs measures in the limit. 

Let Zg_t{x) = E^exp(/? J^v{xs)ds), where is the expectation with respect to the 
measure induced by the Brownian motion starting at x. Thus ^/3,t(0) = Z^^, where Zp^t 
is the partition function introduced in the previous section. 

Theorem 3.1. The process {x{t),0 < t < T} is a time-inhomogeneous Markov process 
with respect to the measures P/3,t- Its transition density is given by 

y)> ^)) = pp(t - s> x)Zp^T-t{x)(Zf)^T-s(y))''^, ■ (5) 
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The transition density q^{{s,y), {t,x)) solves the parabolic equati 

y), {t, x))+ ^Ayq^iis, y), {t, x)) + Vy In Zp,T-s{y)'^yq^{{s, y), {t, x)) = 0. (6) 

With respect to the conditional measure P^,t(- \x{T) — 0), the process {x{t),0 < t < T} 
is a time-inhomogeneous Markov process with transition density 

<lf'^\is,y)dt^x)) = pf^{t- s,y,x)p(^{T -t,x,O){p0{T - s,y,0))-\ (7) 

While this result is not used directly in later sections, it provides some intuition on 
the nature of the limiting processes when we consider the limit T — > oo. 

Proof. The Feynman-Kac formula gives that for < i < T, 

, , s s P3(t,0, x)E^ cxn(3 fl^ ^ v(xs)ds) , , 

P^,T(x(i) e dx) = ' ^ ^^^dx. (8) 

Similarly, for = to < < ^2 < ■■■ < < 2^ and xq — 0, 

P/3,r {x{t}) e dxi, ...,x{tn) e dXn) = 
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117=0 Pf^i^i+i - ti, Xi, Xi+i)E^" exp(/3 /(f *" v{xs)ds) 
— cLx 1(1X2 ■ ■ ■ dXfi . 

Zl3,T 

So, if we set for < s < t < T, 

Qp{{s,y), {t,x)) =Pf}{t - S,y,x)Zf3^T-t{x){Zf),T-s{y))~^, 

then 

n-l 

P^r {x{tl) e dXi, ...,x{tn) e dXn) ^ JJ Qpiiti, Xi), Xj+i)). 

i=0 

Since q^{{s,y), {t,x)) > and 

<ll{{s,y),{t,x))dx^ 1, 

this means that {a;(t),0 < t < T} under the measure Pjs^r is a time-inhomogeneous 
Markov process with transition probabihties . Turning the equation for around and 
solving for yields 

,^ . (lp{{s,y),{t,x))Zp,T-s{y) 

pp{t - s, y, X) = -— . 

Zp,T-t[x) 

Using the fact that 

d 1 

— - s, y, x) + -Aypp{t - s, y, x) + Pv{y)pp{t - s, y, x) = 0, 
we derive that satisfies the equation 

|-«?((^. y), (t, y), (t, x A^'-^-ff 

OS Zp^T-t{x) Zp^T-t{x) 

+ y), (t, x))^^+(3v{yms, yl {t, ^))|^^ 

^ ^P,T-t{X) Zjjj^T-tyX) (9j 

+ r~\ \^l3,T-s{y)+Vyqf3 {{s, y), [t, x)) " 

^ Zp,T-t[X) Zi3^T-t[X) 

=0. 



Simphfying this leads to the following parabolic equation for q^, 

d_ 
ds 



^q^{{s, y), (t, x))+ iA^gJ((s, y), {t, x)) + Vy In Vj,gJ((s, y), (t, x)) = 0. (10) 
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Next we consider the pinned case, for — to < h < ... < tn < tn+i — T and 
Xq — — 0. Then, 

^l3,T{x{ti) e dXi, ...,x{tn) e dXn,x{T) = 0) 



P/3,T {x{ti) e dXi, x{tn) e dXn\x{T) = 0) 



Pf,,T ix{T) = 0) 

Ui=0 Pl3{ti+1 -ti,Xi, Xj+i) 



(11) 



Now set for < s < i < T, 

gf ((s, y), (t, x)) = pfsit - s, y, x)pp{T - t, x, 0){pp{T - s, y, 0))-\ (12) 

Then 



n-l 



P/3,T(a;(ti) e dxi, ...,a:;(tn) e (ia;n|a;(T) = 0) = JJ x^), (t^+i, x^+i)). (13) 

i=0 

Since q^^'^\{s,y), {t,x)) > and 

/ gj((s,y),(^,a;))dx= 1, 

this means that {x{t),0 < t < T} under the conditional measure P^,T(-|a;(T) = 0) is a 
time-inhomogeneous Markov process with transition densities ■ □ 

We shall see below in that in the globular phase f3 > f3cr the drift term V^; In Zp^T-s{x) 
has a non-trivial limit as T oo. This means that for j3 > (3cn the Gibbs measure 
corresponds to a stationary Markov process in the T — > oo limit. On the other hand, 
this hmit will vanish for P < Per- This explains the nature of the diffusive state for high 
temperature. 



4 Critical Value of the Coupling Constant 

Let 

H(i = -A + Pv: L\R'^) L^R'^), v = v{x) e C^{R'^), p > 0. 

We shall always assume that v{x) is non-negative and compactly supported, although 
many results do not require these restrictions or can be modified to be valid without 
these restrictions. We shall also assume that v is not identically equal to zero. It is 
well-known that the spectrum of Hp consists of the absolutely continuous part (— oo, 0] 
and at most a finite number of non-negative eigenvalues: 

aiHp) = (-00, 0] U {A,}, < J < AT, A, = A,(/3) > 0. 

We enumerate the eigenvalues in the decreasing order. Thus, if {Xj} ^ 0, then Aq = 
max Ao. 
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Lemma 4.1. There exists j3cT > (which will he called the critical value of (3) such that 
sup cr(i7^) = for (3 < Per o^nd supcr(if/3) = Ao(/?) > for [3 > Per- For [3 > Per 
the eigenvalue Xo{P) is a strictly increasing and continuous function of p. Moreover, 
Iimi3ii3^^ X{P) =0 and lim/3|oo HP) = oo. 

Proof. The form {Hpip,ip) is positive on a function ip supported on supp(f) if P is large 
enough. Thus sup cr(if^) > for sufficiently large p. On the other hand, cr{Hp) = (— oo, 0] 
when P = 0. Let Per = sup{/3 : sup a{Hf^) = 0}. It is clear that sup a{H/^) = for /3 < Per 
since the operator H/^ depends monotonically on p. 

Other statements easily follow from the fact that for each ip the form {Hpif), ip) de- 
pends continuously and monotonically on p. □ 

Remark. As will be shown below. Per = for d = 1,2, and Per > for d > 3. Thus we 
do not talk about phase transition for ci = 1, 2 since we do not consider negative values 



This implies that there are no eigenvalues for sufficiently small values of /? if (i > 3, that is 
Per > 0. It is also well-know (see [6]) that sup a^Hi^) > for (i = 1, 2 if /? > 0, f > and 
V is not identically zero. These statements will also be proved below without referring to 
the Cwikel-Lieb-Rozenblum estimate. 

5 Analytic Properties of the Resolvent 

The resolvent of the operator Hj^ will be considered in the spaces of square-integrable 
and continuous functions. The resolvent R/siX) = {Hp — A)"-*^ : L^(R'^) L'^iW^) is a 
meromorphic operator valued function on C = C\(— oo,0]. Denote the kernel of Rfi{X) 
by Ri3{X,x,y). If P = 0, the kernel depends on the difference x — y and will be denoted 
by Ro{X,x — y). The kernel Rq{X,x) can be expressed through the Hankel function H^^: 



oip. 



For c? > 3, by the Cwikel-Lieb-Rozenblum estimate [6], 




-Ro(l, x) = c\x\ 2 



lH'^p(zV2\x\), 



(14) 



2 



and 



Ro{X, x) = ck'^-^{k\x\f-iHf {iy/2k\x\), k = y/X, Rek > 0. 



(15) 
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In particular. 



Ro{X,x) 



g-V2fc|x| 



d 



1; Ro{X,x) 



d = 3. 



-V2k ' 



—27i\x 
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We shall say that / G Lg^p(R ) if / is measurable and 



( [ /2(x)el"l'(ix)5 < oo. 



Similarly, we shall say that / G Ccxp{M.'^) if / is continuous and 

lce.p(M<^) = sup(|/(a;)|el"l') < oo. 



Note that i?o(A), A G C, is a bounded operator not only in L^(M'^) but also from 
Cexp(K'^) to C(M'^), where C(]R'^) is the space of bounded continuous functions on W^. 
Denote 

A{\)=v{x)Ro{\): Ll^iR'') Ll^iR") (and Cexp(M'^) ^ C^xp (16) 



The well-known properties of the Hankel functions together with flMj) and f[T5|) imply 
the following lemma (see [8] for a similar statement for general elliptic operators). 

Lemma 5.1. Consider the operator A{X) in the spaces Ll^^iW^) and Cc^piW^). 

(1) The operator A{X) is analytic in X E C . It admits an analytic extension as an 
entire function of VX if d is odd, except d = 1, when it has a pole (with respect to ^/X) at 
the origin. The operator A{X) has the form A{X) = Ai{X) +lnA742(A) if d is even, where 
Ai and A2 are entire functions. 

(2) ^2(0) = if d > 4 (d is even), and therefore A{0) = limA^o.Aec ^(A) exists and 
is a bounded operator for all d> 3. 

(3) The operator A{X) is compact for all A G C U {0} (X^Qifd=l or 2). 

(4) For each e > 0, we have ||A(A)|| = 0(1/|A|) as X 00, |argA| < n — e. 

(5 ) The operator A{X) has the following asymptotic behavior as A — > 0, X E : 

A{X) = -vPi/^/X + 0{l), rf=l, 
A(A) = -t;P2ln(l/A) + 0(l), = 2, 
A{X) = -v{Ps + Qs^/X) + 0(|A|), d = 3, 
A{X) = -v{Pi + Q^X ln(l/A)) + 0(|A|), = 4, 
A{X) = -v{Pd + QA) + 0(\Xf^), d>5, 
where the operators Pd, d > 1, Qd, d > 3, have the following kernels: 

^i(^) y) = ^2(a;, y) = -, 

^3(2;, y) = 7—^ r, Qsix, y) = — 

27r\x-y\ V27r 

1 1 

P4{x,y) = — , Q4{x,y) 



= \x-l\<i-^ ' = ^d~4)\x-y\''-^ ' ad>0, d>5. 
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Proof. Let d be odd. From f|T^ . (fT5|) and (fT6l) it follows that the kernel A{X,x,y) = 
v{x)Ro{X, X — y) oi the operator A(X) is an entire function of A; = ii d > 3 (but has a 
pole at = if d = 1). The kernel has a weak singularity aX x = y and an exponential 
estimate at infinity. To be more exact, 

\A{e, x,y)\ + l^^i^l^l < C{d, kMx)\e^''^^-y^\\x - y\ + \x - yr^'^-^)), (17) 
were C{d, k) has a singularity at = if ci = 1. Since 

P(A:^)||c..p(M^) < sup / e\^\'~\y\'\Aie,x,y)\dy, 



< sup / el^l "1^1 I ''' \dy. 



the estimate f|T7|) immediately leads to the analyticity in A; = -\/A of the operator A{X) in 
the space Cexp(IR'^)- In order to get the same result in the space Lg^p(R'^), we represent 
A(X) in the form Bi + B2 were the kernel Bi{X,x,y) of the operator Bi is equal to 
x{x — y)A{X,x,y). Here x is the indicator function of the unit ball. Since 

\x{x)Ro{k',x)\ + \^x{x)Mk',x)\ G L\R^), (18) 

the convolution with x{x)Ro{k'^, x) is an analytic in k operator in the space L^(M'^). Then 
Bi (which is the convolution followed by multiplication by v{x)) is an analytic operator 
in the space Lg^p(M'^). The product of the kernel of the operator B2 and e'^^l "'^1 is square 
integrable in {x, y). The same is true for the derivative in k of the kernel of B2 multiplied 
by gki -\y\ Thus B2 is also analytic in k. This completes the proof of the analyticity of 
74(A) when d is odd. The case of even d is similar. One needs only to take into account that 
Ro{X,x) has a logarithmic branching point at A = in this case. The second statement 
of the lemma follows immediately from f|T4|) . f|T5l) and f|T6l) . 

To prove the compactness of A{X), we note that the estimate f|T7j) is valid not only 
for A{k'^, X, y) and dA{k'^, x, y)/dk, but also for VxA{k'^, x, y). Thus the arguments above 
lead to the boundedness of the operators -^A^X) (the composition of A{X) with the 
differentiation). Since the supports of functions A{X)f belong to the support of f, the 
standard Sobolev embedding theorems imply the compactness of the operator A{X) in 
both the spaces Ll^^iM.'^) and Cexp(M^). 

In order to prove the fourth statement of the lemma, we observe that the L?{W^) norm 
of the resolvent Ro{X) does not exceed l/|ImA| (the inverse distance from the spectrum). 
Since A{X) is obtained from -Ro('^) after multiplying it by a bounded function with compact 
support, the Lg^p(R°') norm of A{X) does not exceed c/|ImA|, where c is a positive constant 
which depends on v. The norm of A{X) in the space Cexp(K'^) can be estimated by 
sup^gjgd |f (x)el^l I Jjgd \Ro{X,x)\dx., which is of order 0(1/|A|) as A — 00, |argA| < vr — £, 
due to ([15]). 
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The remaining statements also easily follow from ([Hj) and (fTSjl . □ 
Note that for d > 3, there exists the limit 

i?o(0, x-y) := ^Jim ^, ^o(A, x - y) = -a^lx - y\'^~'^, 

which is a fundamental solution of the operator |A. The operator with this kernel will 
be denoted by Ro{0). While Ro{X), A G C, acts in L^iR'^) and C(M'^), the operator Ro{0) 
only maps Cexp(IR'^) to C(M'^) if (i < 5. The following lemma follows from formulas (HM 
and (fT5|) similarly to Lemma 15.11 

Lemma 5.2. Ford > 3, the operator Rq{X) considered as an operator from Cexp(R°') to 
C(M"') is analytic in X & C It is uniformly bounded in C. For each e > 0, it is of order 
0(1/|A|) as X — * oo, |argA| < n — e. It has the following asymptotic behavior as X 0, 
X e C: 

i?o(A) = i?o(0) + 0(VW), d = 3, 

i?o(A) = i?o(0) + 0(|AlnA|), d = 4, 

Ro{X) = Ro{0) + 0{\X\), d>5. 

The following lemma is simply a resolvent identity. It plays an important role in our 
future analysis. 

Lemma 5.3. For X G C, we have the following relation between the meromorphic 
operator-valued functions 

Rp{X) = i?o(A) - i?o(A)(/ + (3v{x)Ro{X))-'[Pv{x)Ro{X)] (19) 

Remark. Note that (fT9|) can be written as 

R(,{X) = RoiX) - i?o(A)(/ + PAiX))-'[Pv{x)RoiX)]. (20) 

From here it also follows that 

R(s{X) = RoiX){I + pA{X))-\ (21) 

which should be understood as an identity between meromorphic in A operators acting 
from L2^p(M'^) to L'^{R'^) and from Cexp(M^) to C(M^). In the lattice case considered in [1], 
the operator A{X) has rank one and 

Rf,{X,x,y) = Ro{X,x,y)/{l - PI{X)), 

where /(A) is an analytic function of related to ^(A). This exact formula is the key 
to all the results in |T]. 

The kernels of the operators / + l3A{X) (both in spaces Ll^^iM.'^) and Ccxp(IR'^)) are 
described by the following lemma. 
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Lemma 5.4. (1) The operator-valued function {I + (]A{X))^^ is meromorphic in C. It 
has a pole at X & C if and only if X is an eigenvalue of Hp. These poles are of the first 
order. 

(2) Let Xi{P) be a positive eigenvalue of Hp. There is a one-to-one correspondence 
between the kernel of the operator I + i3A{Xi) and the eigenspace of the operator Hp 
corresponding to the eigenvalue Aj. Namely, if (J + (3A{Xi))h = 0, then ip = —Ro{Xi)h is 
an eigenfunction of Hp and h = jSvip. 

(3) If d > 3, there is a one-to-one correspondence between the kernel of the operator 
I + (3A{0) and solution space of the problem 

1 dib 

Hpm = -Aij + Pv(xU = 0, ^(x) = Oflxp-'^), -^(x) = 0(\x\^-'^) as r=\x\^oo. 

2 or 

(22) 

Namely, if (/ + (3A{Q))h = for h E L^^piM'^), then h e Cexp(ffi'^); ^ = -i?o(0)/i is a 
solution of and h = Pvip- 

Remark. The relations fl22|) are an analogue of the eigenvalue problem for zero eigenvalue 
and the eigenfunction which does not necessarily belong to L'^{W^) (see Lemma 15.61 
below). We shall call a non-zero solution of (1221) a ground state. 

Proof. The operator A{X), X G C, is analytic, compact, and tends to zero as A ^ +oo by 
Lemma IHTTl Therefore {I + f3A{X))~^ is meromorphic by the Analytic Fredholm Theorem. 

If A e C is a pole of (/ + PA{X))~^, then it is also a pole of the same order of Rp{X) 
as follows from (12T!) since the kernel of -Ro(A) is trivial. Therefore the pole is simple and 
coincides with one of the eigenvalues Aj. Note that A is a pole of {I + j3A{X))^^ if and only 
if the kernel of / + PA{X) is non-trivial. Let h G L'^^^iM.'^) be such that ||/?.||l2^ (Rd) 7^ 
and (/ + (3vRo{X))h = 0. Then ip := -Ro{X)h G ^^(M^) and (^A - A + Pv)ij = 0, that is 
ip is an eigenfunction of Hp. 

Conversely, let ip G L^(M'^) be an eigenfunction corresponding to an eigenvalue Aj, 
that is 

{^A-X,)ij + Pvi: = 0. (23) 

Denote h = j3vip. Then (^A — Xijip = —h. Thus = —Ro{Xi)h and fl23l) implies that h 
satisfies (/ + PvRo{Xi))h = 0. Note that h G C°°(M'^), h vanishes outside supp(v), and 
therefore belongs to the kernel of / -|- pA{Xi). This completes the proof of the first two 
statements. 

Similar arguments can be used to prove the last statement. If /i G L'^^p{M.'^) is such 
that ||/i||L2^ (iRd) 7^ and (J -|- f3A{0))h = 0, then h has compact support and the integral 
operator -Ro(O) can be applied to h. It is clear that ip := —RQ{0)h satisfies (!22l) and, since 
h has compact support, h G Cexp(K'^)- 

In order to prove that any solution of (1221) corresponds to an eigenvector of I + l3A{0), 
one only needs to show that the solution of the problem fl22|) can be represented in the 
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form ip = — i?o(0)/i with h = f3vip. The latter follows from the Green formula 

ilj{x) = -(i?o(0)(/5t;^))(a;) + / [i?o(0,x - |-i?o(0, x - |x| < a, 

J\y\=a or 

after passing to the limit as a — > cxd. □ 
Lemma [5.41 can be improved for A = Ao(/3). Due to the monotonicity and continuity 
of A = Ao(/3) for [3 > Per, we can define the inverse function 

/3 = /3(A):[0,oo)->[/?,„oo). (24) 

We shall prove that the operator —A{X), A > 0, has a non-negative kernel and has 
a positive simple eigenvalue such that all the other eigenvalues are smaller in absolute 
value. Such an eigenvalue is called the principal eigenvalue. 

Lemma 5.5. The operator —A{X), A > 0, (in the spaces Lg^p(R'^) and CcxpiM'^)) has the 
principal eigenvalue. This eigenvalue is equal to 1//3(A) and the corresponding eigenfunc- 
tion can be taken to be positive in the interior of supp(f ) and equal to zero outside of 
supp(f). If d>3, then the same is true for the operator —A{0) (in particular, f3cr > 0^). 

Remark 1. Let d > 3. Lemmas 15.41 and 15.51 imply that the ground state of the operator 
Hp for [3 = Per (defined by fl22l) ) is defined uniquely up to a multiplicative constant and 
corresponds to the principal eigenvalue of A{0). The ground state (with A = 0) does not 
exist if P < Per- 

Remark 2. Let d>3. From Lemma [5.11 it follows that 

lim A(X) = A(0). 

a^o.agc 

Therefore for all A G C with |A| sufficiently small, the operator —A{X) has a simple eigen- 
value whose real part is larger than the absolute values of the other eigenvalues. We shall 
denote this eigenvalue by 1//3(A), thus extending the domain of the function /3(A) (see 
flM|) ) from [0, oo) to [0, oo) U (f/flC), where [/ is a sufficiently small neighborhood of zero. 

Proof of Lemma \5.5\ By Lemma 15.41 it is sufficient to consider the case of Lg^p(]R'^). 
The maximum principle for the operator (^A — A), A > 0, implies that the kernel of the 
operator Ro{X), A > 0, is negative. Thus, by (fT6|l . for all y the kernel of —A{X) is positive 
when X is in the interior of supp(w) and zero otherwise. Thus —A{X), A > 0, has the 
principal eigenvalue (see [4]). On the other hand, by Lemma 15.41 1/P{X) is a positive 
eigenvalue of —A{X). Note that this is the largest positive eigenvalue of —A{X). Indeed, 
if fi = 1/P' > 1/P{X) is an eigenvalue of ~A{X), then A is one of the eigenvalues Aj of Hp' 
by Lemma [5.41 Therefore, Xi{P') = Xq{P) for P' < p. This contradicts the monotonicity 
of Xq{P). Hence the statement of the lemma concerning the case A > holds. 

For d > 3, the kernel of —^(0) is equal to vPd and has the same properties as the 
kernel of —A{X), A > 0. Thus —^(0) has the principal eigenvalue. Since A{X) — > A{0) 
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as A I 0, the principal eigenvalue 1//3(A) converges to the principal eigenvalue /i < oo of 
—A{0). On the other hand, /5(A) is a continuous function, and therefore /x = l/Pcr, which 
proves the statement concerning the case A = 0. □ 

The relationship between ground states and eigenfunctions of Hp is explained by the 
following lemma. 

Lemma 5.6. Let (3 = Per- If d = 3 or d = 4, then Hp has a unique ground state (up to a 
multiplicative constant), hut \ = Q is not an eigenvalue. If d> 5, then X = is a simple 
eigenvalue of Hp and the sets of ground states and eigenfunctions coincide. 

Proof. The ground states belong to L^(M'^) if and only ii d > 5. In order to complete the 
proof of the lemma, it remains to show that any eigenfunction of Hp with zero eigenvalue 
satisfies ( 122|) . Thus, it is enough to prove that if ^^ip + Pv{x)'iIj = and ^ L'^iM.'^), 
then ip = — i?o(0)/i with h = j3vip. From ^A?/' + f3v{x)ip — Xip = —Xip we obtain ip = 
-Ro{X){h + Xip). Obviously RoiX)h RoiO)h in L'^{R'^) as A j since h G LLpl^'^)- 
Now the lemma will be proved if we show that 

l|Ai?o(A)^||i.(«.) = jj^^^fdo - as A i 0. 

The latter follows from the dominated convergence theorem. □ 
The following lemma summarizes some facts about the operator [I + (3A{X))~^ proved 
above. It also describes the structure of the singularity of the operator (J + [3A{X))~^ for 
A and /3 in a neighborhood of A = 0, /3 = /5cr. 

Lemma 5.7. Let > 3 and /3 > 0. The operator (/ + I3A{X))~^ (considered in Lg^p(R'^) 
and Cexp(K'^) ) is meromorphic in X & C and has poles of the first order at eigenvalues of 
the operator Hp. For each e > and some A = A{(3), the operator is uniformly bounded 
in X e C, |argA| < n — e, \X\ > A. 

If f3 = Per, then the operator {I+f3A{X))~^ is analytic in X & C and uniformly bounded 
in X e C, |argA| < n — e, \X\ > e. 

If (3 < Per, then the operator {I+PA{X))^^ is analytic in X E C and uniformly bounded 
in X E C , |argA| < n — e. 

There are Aq > and 6o > such that for A G C U {0}, |A| < Aq, \P — Pcr\ < ^O; 
f3 7^ P{X), we have the representation 

(J + /3AiX))-' = fj^^iB + S,{X)) + C(A, P). (25) 

Here /5(A) is defined in Remark 2 following Lemma l575\ B is the one dimensional operator 
with the kernel 
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where ip is a ground state defined in the Remark following Lemma 



and 



Ssi\)=0{^/\\\), 54(A) =0(|Aln(A)|), Sdi\)=0{\\\), d>5, as A ^ 0, A e C, 

(27) 

Sd{0) = 0, d > 3, and C(A,/3) is bounded uniformly in A and [3. 



Proof. The analytic properties of (/ + I3A{X))^^ follow from Lemma [5.41 By Lemma [5. 11 
the norm of A{X) decays at infinity when A 00, |argA| < vr — £. Therefore there is 
A > such that the operator (/ + [3A{\))~^ is bounded for |argA| < vr — |A| > A. 

If /5 < (3cr, then (/ + (3A{\))~^ does not have poles in A G C, and therefore A can be 
taken to be arbitrarily small. 

\i (3 < Per, then {I + i3A{0)) is invertible by Lemma [531 By Lemma [STTl the operators 
A{X) tend to A{0) when A ^ 0, A G C. Therefore {I + (3A{X)y\ A G C, are bounded in 
a neighborhood of zero. It remains to justify fl25l) . 

For d > 3, let hx be an eigenvector corresponding to the eigenvalue 1//3(A) of the 
operator —A{X), A G [0, 00) U {U (1 C). By Lemma [5.51 and the second remark following 
it, this eigenvector is defined up to a multiplicative constant. Let A*{X) be the opera- 
tor in Lexp(^'^) or C'exp(K'^) with the kernel A*{X,x,y) = A(A, x)el2^l'-l^l'. Similarly to 
Lemma [5751 it is not difficult to show that 1/(3{X) is an eigenvalue for the operator —A*{X) 
and that its real part exceeds the absolute values of the other eigenvalues. The corre- 
sponding eigenvector h*^ is uniquely defined up to a multiplicative constant. Moreover, 
we can take hx and h\ such that 

t;(x)el^l'/i*(x) = hxix). (28) 

Note that hx and h\ can be chosen in such a way that 

\\hx-ho\\,\\h*x-hl\\<k\\AiX)-A{0)\\ (29) 

for some k > and all sufficiently small |A|, where the norms on both sides of (!29l) are 
either in the space Ll^pO^'^) or Ce^p{R'^). 

Recall that A{X) —>■ A{0) as A — > 0, A G C, by Lemma [STTl Using this and the fact that 
1/Pcr is the principal eigenvalue for —A{0), it is easy to show that there are Ai > and 
5i > such that for A G C U {0}, |A| < Ai, the eigenvalue 1//3(A) of the operator — A(A) 
is the unique eigenvalue whose distance from I /Per does not exceed 61. Take < Aq < Ai 
and < 60 < 61 such that for A G C U {0}, |A| < Aq, the distance between 1//3(A) and 
I /Per does not exceed 60. 

Then for A G C U {0}, |A| < Aq and P such that \1/P — l/Per\ < Sq, the operator 
valued function 



is meromorphic inside the circle 7 = {z : \z—l/Per\ = It has two poles: one at z = 1/P 
and the other at z = 1/P{X). The residue at the first pole is equal to {A{X) + I/P)^^. In 
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order to find the residue at the second pole, recall that it is a simple pole for {A{\) + zI) ^, 
and therefore 

(A(A) + zir' = r_i(A)(^ - + ro(A) + Ti(A)(^ - ^) + - 

for some operators T_i,Tq,Ti, ... and all z in a neighborhood of 1//3(A). From here and 
the fact that the kernels of A{X) +///?( A) and A*{\) + J//3(A) are one-dimensional and 
coincide with span{h\} and span{/i^}, respectively, it easily follows that 



r_i(A)/= feLl^iR") (in particular if /Ga.p(M'^)). 



From ([29D and Lemma O it follows that Sd{\) := T_i(A) - T_i(0) satisfies The 
residue of F{z) at z = 1//3(A) is equal to 

^(T..(0) + ..(A)), 

Integrating F{z) over the contour 7, we obtain 

The right hand side of this formula is uniformly bounded, which completes the proof of 
the lemma if we show that T_i(0) = B. Thus it remains to prove that 

ho{x)e\y\"hl{y) _ v{x)tlj{x)^{y) 

The latter follows from the relation ho = f3vip (see Lemma [5.41) and ( l28l) . □ 
Formula (12T|) and Lemmas 15.21 and 15.71 imply the following result . 

Lemma 5.8. Let d > 3 and (3 > 0. The operator RpiX) (considered as an operator 
from Cexp(IR'^) to C{M.'^)) is meromorphic in X E C and has poles of the first order at 
eigenvalues of the operator Hp. For each £ > and some A = A(/3), the operator is 
uniformly hounded in X E C , |argA| < n — s, |A| > A. It is of order 0(1/|A|) as X ^ 00, 
|argA| < 71 — e. 

If P = Per, then the operator Rfs^X) is analytic in X E C and uniformly bounded in 
X E C, |argA| <7r-e, |A| > e. 

If (3 < Per, then the operator RfsiX) is analytic in X E C and uniformly bounded in 
X E C, |argA| <TT-e. 

There are Aq > and 60 > such that for X E C , < \X\ < Xq, \P — Pcr\ < ^O; 
P ^ P{X), we have the representation 

^/^(^) = + ^'^(^)) + (30) 

where /5(A) is defined in Remark 2 following Lemma l575\ and B is given by ^2B^) . Sd, d >3, 
satisfy [21^ , and C{X,P) is bounded uniformly in X and p. 
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6 The Behavior of the Principal Eigenvalue for P | per 

In Lemma 15.51 we showed that (3cr > for d > 3. The following theorem implies, in 
particular, that (3cr = for li = 1 or 2. 

Theorem 6.1. For d = 1,2 (when (3cr = 0^ the eigenvalue \o{(3) has the following 
behavior as [3 [ Per- 

^oiP) ~ -CiP'^, ci = I v{x)dx, d = 1, (31) 

Ao(/?)~exp(-|), C2 = -, d = 2. (32) 
In dimensions d > 3 the eigenvalue Xo{(3) has the following behavior as (3 I (3cr- 

Ao(/5) ~C3(/?-/5cr)', d = 3, (33) 

Ao(/?) ~ C4(/5 - /?,,)/ ln(l/(/3 - /?,,)), = 4, (34) 
Ao(/3)~q(/3-/5,,), d>5, (35) 
where q 7^ 0, ci > 3, depend on v and will be indicated in the proof. 

Proof. Since we are interested in the behavior of \o{P) for P | P^r and \o{P) I when 
/? i Per by Lemma I^TTl we shall study the behavior of /?(A) as A | (or, more generally, 
as A — > 0, A G C). The arguments below are based on Lemma [5. 1[ 

First consider the case d = 1. For A — > 0, A G C, the eigenvalue problem for —A{X) 
can be written in the form 

{vPi + 0{^/X))hx = ^hx. (36) 

Note that the kernel of vPi is positive when x is an interior point of supp(f). Therefore 
vPi has a principal eigenvalue. In fact, the operator vPi is one-dimensional and the 
eigenvalue is equal to Ci/a/2 where Ci = J^av{x)dx. Since this eigenvalue is simple and 
the operator in the left-hand side of (|36ll is analytic in V^, both hx and VX/ P{X) are 
analytic functions of \/A in a neighborhood of the origin and 

lim (y/X/p(X)) = C1/V2. 

Therefore, Per = 0, P{X) is analytic in v^A, and P{X) ~ v^2A/ci as A — 0, A G C, which 
proves (13T|) . 

The same arguments in the case d = 2 lead to the relation 

lim ( , I , ) = ci/n. 
A^o,AeC'7?(A) InA 
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This implies that Per = and (l32l) holds. 

In the case d = 3 the eigenvalue problem for —A{X) takes the form 

(-A(O) + VXv{x)Q3 + 0(A))/lA = -r^/iA- (37) 

As in the one-dimensional case, 1//?(A) and hx are analytic functions of V^- Now l/(3cr is 
equal to the principal eigenvalue of ~A{0). Recall that is the principal eigenfunction of 
—A{0) and /iq is the principal eigenfunction of —A*{0). Standard perturbation arguments 
imply that 

J_ = i__,yX + 0(A). A^O.A.C, (38) 

where 

^^-W, (39) 

which implies (ESD with cg = 1/(72/?^^). Note that 7 > since the kernel of the operator 
VQ3 is negative and principal eigenfunctions ho, can be chosen to be positive inside 
supp(w). 

Formula for 7 can be simplified. We choose ho = jSvip (see Lemma [5.41) and defined 
in ([28]). Then 

7 = V ;f-v ; ; ^ ^ _ 3_ ^^^^ 

V27r jjg3 f (x)V'^(x)dx 

Let (i = 4. Then instead of flTTI) we get 

(-A(O) + A ln(l/A)t;Q4 + 0{\))hx = -^h^. (41) 
From here it follows that 

^ - 7A ln(l/A) + 0(A), A 0, A G C, (42) 



/3(A) 

where l//3cr is the principal eigenvalue of —^4(0) and 7 is given by fl5I?]) with Q3 replaced 
by Q4. Thus ^ holds with C4 = l/inl^lr). 
For (i > 5 we get 

(-A(O) + At;Qd + 0(A3/2))^^ ^ 1 
From here it follows that 

_l_._l_,, + 0(An A^O.A.C, 

where l//3cr is the principal eigenvalue of —^4(0) and 7 is given by fl39|) with Qs replaced 
by Qd. Thus ^ holds with q = 1/(7/?^^,). □ 
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7 Asymptotics of the Partition Function, Solutions, 
and Fundamental Solutions 

We shall need the following notation. Recall from (jlj) that by pfs{t,y,x) we denote the 
fundamental solution of the parabolic problem 

g-^ = ^^xPisit, y, X) + f3v{x)pf^{t, y, x), 

pp{0,y,x) = 6{x-y). 

For a given / G ^^(R'^), let 

up{t,x)= pp{t,y,x)f{y)dy 

be the solution of the Cauchy problem with the initial data /. The partition function is 
defined as the integral of the fundamental solution 

Zfs^x) = / p,3{t,x,y)dy = / pp(t,y,x)dy. 

Note that the partition function defined in ([3]) is simply Zf3^T = ^/3,t(0). Also note that 
Zf^^x) is the solution of the Cauchy problem with initial data equal to one: 

= ^AZp,{x) + (3v{x)Zp,{x), Z^,o(x) ^ 1. 

For P > Per, let V'/? be the positive eigenfunction for the operator Hp with eigenvalue 
Ao(/5) normalized by the condition | l^/'/j] |l2(r) = 1. This function is defined uniquely by 
Lemma and is equal to —Ro{X)hx, where A = Ao(/9) and hx is the principal eigenfunc- 
tion for the operator —A{X). Note that ipf3 decays exponentially at infinity. 

For a G M, let T{a) be the following contour in the complex plane 

r(a) = {a — s + is, s > 0} U {a — s — is, s >0}. 

We choose the direction along T{a) in such a way that the imaginary coordinate increases. 
The following lemma is an important tool for investigating the asymptotics of ^/3,r- 

Lemma 7.1. Let a > Ao(/3). Then for f G L'^{R'^) (or f G C^^piR"^) ) and t > 0, 

up{t,x) = ^.f e^\Rp{X)f){x)dX, (43) 

which holds in L^(R'^) (or C^W^)). This formula remains valid if the initial function f 
is identically equal to one and R[^{X)f is understood by substituting f = I into ^^) with 
i?o(A)l = — 1/A. More precisely, 

-,\t 



Zp,{x) -1 = 7^ [ ^—{Rp{X){(3v)){x)dX (44) 
^7r« Jr{a) X 

m L\R'^) and C{R'^). 
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Proof. First, let / G L^(M'^). We solve the Cauchy problem for Ujs using the Laplace 
transform with respect to t. This leads to (143!) with r(a) replaced by the line {A : ReA = 
a}. The integral over this line is equal to the integral over T{a) since the resolvent is 
analytic between these contours and its norm decays as |A|~^ when |A| — > oo. 
Now let / = 1. Then w{t, x) = Zp^t{,x) — 1 is the solution of the problem 

^^1^' = -Au;(t, x) + (3v{x)w{t, x) + (3v{x), w(0, x) = 0. 
By the Duhamel formula and fH3l) . 



w(t,x) = — I I e^^'-'\R3{X)pv)ix)d\ds 

27n Jn ./r(a) 



/ '-^{R,{X)Pv){x)dX = / ^-!^{R,{X)Pv){x)dX, 
27rz Jr^a) X Im J^^^^ A 

since in the domain T^{a) to the right of the contour T{a), the operator -R/3(A) : L^(]R^) 
L^(R'^) is analytic and decays as lAI^""^ at infinity. This justifies (jH]) in L^(M'^) sense. It 
remains to show that the right-hand side of fH51) is continuous for / G Ccxp(K'^) and the 
right-hand side of (jH]) is continuous. Since l3v G C^, the integrands are continuous in 
{t,x) for each A G F^. It remains to note that the integrals converge uniformly when 
X G M", t>to>0. This is due to the fact that \ \Ri3{X)f\\c(Rd), \ \R[^{X)Pv\\c(Rd) < Cd{a), 
as follows from Lemma [5.81 □ 
In order to state the next theorem we shall need the following notation. As in part 
(3) of Lemma [531 it is not difficult to show that for > 3, < /3 < per and / G C^Ir'^) 
there is a unique solution of the problem 

Hi3{^) = ^A^ + Pv{x)ip = f, (^ = 0(|x|2-^), ^(x) = Oi\x\'-'') as r = \x\ ^ oo. 

(45) 

This solution is given by = RoiO){I + (3A{0)y^f. For / = -(3v, we denote this solution 

by <^/3. 

Theorem 7.2. (1) For j3 > /Scr there is e > such that we have the following asymptotics 
for the partition function: 



Zf3,ti^) - 1 = exp(Ao(/3)t)(||^/3||ii(Kd)?/'/3(x) +0(exp(-£:t))) as t 



oo, 



which holds in L'^{W^) and in C(M'^), where ipp is the positive eigenf unction for the operator 
Hp with eigenvalue Ao(/3) normalized by the condition | IV'/?! |l2(k) = 1. 

(2) For 13 = Per we have the following asymptotics for the partition function: 

Zf3,t{x) = k^t^^'^i/j{x) + 0(1) as t ^ oo, d = 3, 

Zf3,ti^) = f^^T-i^i^) + OirA-) as t OO, d = A, 
mr in t 
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Zi3,t{x) = kdtilj{x) + 0{Vi) as t ^ oo, d > 5, 

which holds in C(M.'^). Here kd, d> 3, are positive constants and ip is the positive ground 
state for Hp^^ normalized by the condition \ \i3crVip\\L2^ (^d) = 1. 
(3) IfO< (3 < Per, then 



lim Zfs^tix) = 1 + (pp{x) 

t— >oo 



in C(l 



Proof. (1) Note that the resolvent RpiX) has only one pole between the contours r(a) and 
r(Ao(/?) —e)iie is less than the distance from Aq to the rest of the spectrum. This pole is 
at the point Ao(/3) and the residue is the integral operator with the kernel —iIj p{x)ip piy) . 
Therefore from (jHj) it follows that 

ZfsA^) - 1 = ^-j^M^) / My)My)dy - tt- / —iRf3iX)i3v)ix)dx. (46) 

AolP) jRd -^vrz Jr(Xo{f3)-s) A 

Since (|A + Pv — \o{P))ipi3 = 0, we have Pvip/s = (Ao(/?) — ^A)ipf^, and the integral in 
the first term of the right-hand side of (H^ is equal to Ao(/5)||V'/3||Li(R'i)- Thus the first 
term on the right-hand side coincides with the main term of the asymptotics stated in 
the theorem. 

It remains to show that the second term on the right-hand side of fj46l) is exponentially 
smaller than the first term. This is due to the fact that the norm of the operator R/^iX) 
is of order 1/|A| at infinity for A G r(Ao(/9) — s). 

(2) Let d = 3. First, let us analyze when p = p^r and A 0, A G C. By 
the factor P{X)/{p{X) - P) in the right hand side of ([30]) is equal to {Pcrl\f\y^ + 0(1) 
as A — 0, A G C, where 7 > is given by (1391) . 

We choose the same ground state %p specified in the statement of Theorem 17. 2[ Then 
from (12^ and Lemma 15.41 it follows that 

, . , X fadV(x)tp(x)dx , ^, , fj,av(x)ih(x)dx , , 

jj^a'v{x)tp'^[x)ax jj^av{x)tp'^[x)ax 

Now, by Lemma EJ] and ([33]), (EHD, 

R^JX)iPerv) = ^ + D(A) = zM + DiX), 4>0, (48) 

'yPcrVX jj^avix)i>'^ix)dx VA 

where the remainder D{X) is of order 0(1) when A — > 0, A G C Note that -D(A) is 
bounded on r"'"(0) since the left hand side and the first term on the right hand side of 
are bounded on r"'"(0) outside a neighborhood of zero. 
Next, we apply ([S]) with a replaced by 1/t and use the expression fllS]) to obtain 



1 r e^* ^ k'r^^ 
27ri Jr{i/t) X y/x 



ZpA^) - 1 = / ^(4^ + D{X))dX. (49) 
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Let us change the variables in the integral Xt = z. Thus 

1 f e%Vik'ip 



2Tii Jy[i) z ^Jz t 

The contribution to the integral from the term containing D{z/t) is bounded, while the 
contribution from the first term is equal to k^t^^'^ip^x), as claimed in the lemma. One 
needs only to note that > since 

— / z~'/'e^dz = - [ z-'/'e^dz = - T a-'/'e-'^da = ^ > 0. 
2vr« Jr(i) TTi yr(i) tt a/tt 

If = 4, then ([MD, imply that /5(A) - Per ~ P^-fXHl/X) as A ^ 0, A G C This 
leads to the following analog of (149|) 

Z,A^)-1 = ^[ '^i^f^, + D{X))dX, k',>0, 
27rz Jr(i/t) A Aln(l/A) 

where -D(A) is of order 0(l/|Aln^A|) when A — > 0, A G C and is bounded at infinity. 
After the change of variables Xt = z, we obtain 

Zp^tix) - 1 = — / —[^7— — ] — T + D{-))dz, 
2m Jr{i) -2 -2(lnr — mzj t 

which easily leads to the second part of the lemma in the case d = 4. The treatment of 
the case d > 5 is similar. 

(3) We apply with a replaced by 1/t to obtain 



Z^,{x) -1 = 7^ [ ^-^{Rp{X)if3v)){x)dX = / ^-l{R^d){f]v)){x)dz. (50) 



'r(i/t) ^ Jr(i) 
Note that by Lemma [5.21 and since is not an eigenvalue of ^(0) we have 



\im Rf,{X){f3v) = \im Ro{X){I + f3A{X))-\(3v) = i?o(0)(/ + /3A(0))-i(/3t;) = -y.^. 

Since the difference between Rp{z /t){l3v) and —ipjs is bounded on r(l), one can pass to 
the limit t ^ oo under the integral sign in fl50|) . which leads to 



hm Zp,{x) = 1 + ^ / -dz = l + 
t^oo 2m Jr(i) z 



iPf3[X) 

(1) 

□ 

The third part of Theorem 17.21 establishes the existence of limt^oo Zp^tix) for (3 < (3cr- 
Next we examine the behavior of this quantity as (3 ] Per- 
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Lemma 7.3. There are positive constant bd, d > 3, such that 

\imZp^t{x)-l = —^ilj{x)+0{l) as (5 ] f5,r 
is valid in C(M.'^), where ip is the positive ground state for Hp^^ normalized by the condition 

Proof. By the third part of Theorem 17.21 we only need to find the asymptotics as /5 | (3cr 
of (pf3 = -i?o(0)(/ + l3A{0))-\(3v). From ([25]) with A = and /3(0) = per and it 
follows that 

= -i?o(0)(/ + PA{0))-\Pv) = .^±^R,(Q)B{M + 0(1) = + 0(1) 

for some positive constant h^- □ 



8 Behavior of the Polymer for [3 > per 

In this section we shall assume that f3 > (3cr is fixed. A result similar to the first part 
of Theorem 17.21 is valid for the solution of the Cauchy problem and for the fundamental 
solution. 

Theorem 8.1. Let f e L'^{R'^) (or f G Ccxp(M'^)| For (3 > there %s e > such 
that we have the following asymptotics for the solution up of the Cauchy problem with the 
initial data f : 

Ufs{t) = exp(Ao(/5)t)((V^/3, f)L^R'i)i^t3 + qf{t)), (51) 

which holds in L^(M'^) (or in C(]R'^)j, where \ \qf{t)\\ < c| |/| | exp(— et) for some c and all 
sufficiently large t. 

We have the following asymptotics for the fundamental solution of the parabolic equa- 
tion: 

Pl3it,y,x) = exp{Xo{(3)t){ipf3{y)ip^{x) +q{t,y,x)), (52) 

where limj^oo ||q'(''^;1/;^)|| = 0, uniformly in y, and (Ell) holds in L^(M.'^) and in C{R'^) for 
each y fixed. 

Proof. The proof of (ISTll is the same as the proof of the first part of Theorem 17.21 and 
therefore we omit it. 

Let /^'^(x) = pp{5, y,x) be the fundamental solution of the parabohc problem at 

time 6. Note that /^'^ G L'^{R'^) for all (5 > and all y, and f^/ G Cexp(M'^) for all 
sufficiently small S > and all y. Denote the solution of the parabolic equation with the 
initial data /^'^ by M^'^(t,x). Then 

Ppit, y, x) = uf{t -6,x)= exp(Ao(/5)(t - fp^) L^Rd)'^t3{x) + q\t, y, x)), 



23 



where | x) 1 1 < c||/j'^|| exp{—e{t — 6)) for some c and all sufficiently large t. 

Note that {tpis, f^^) L'^(^d-^ can be made arbitrarily close to 'ippiy) uniformly in by 
choosing a sufficiently small 5, and 1 | is uniformly bounded in y for any fixed 6 . This 
justifies (E2D. □ 

Next, let us study the distribution of the end of the polymer with respect to the 
measure P/j^t as T — oo. 

Theorem 8.2. The distribution of x{T) with respect to the measure Pp^r converges, 
weakly, as T ^ oo, to the distribution with the density 'ip/s/W^'/sWLHRd)- 



Proof. The density of x(T) with respect to the Lebesgue measure is equal to 
Z^,t(0) exp{\o{P)T){\\MLHR^)MO) + o{l)y 



(53) 



where q is the same as in (!52|) . When T — oo, the right hand side of (153|) converges to 
■^z?!^)/! I"^/?! IlHik'*) uniformly in x by Theorem 18. 1[ This justifies the weak convergence. □ 
Now let us examine the behavior of the polymer in a region separated both from zero 
and T. Let S{T) be such that 

lim S{T) = lim (T - S{T)) = +oo. (54) 

T — >oo T — >oo 

Let s > be fixed. Consider the process y^(t) = x{S(T) + 1), < t < s. 

Theorem 8.3. The distribution of the process y'^{t) with respect to either of the measures 
Pf3,T or Pj3^T{-\x{T) = 0) converges as T ^ oo, weakly in the space C{[0, s],M.'^) , to the 
distribution of a stationary Markov process with invariant density ip"^ and the generator 

2 " i'g 

Remark. Let 

rpit, y, x) = ^^-^^^ exp(-Ao(/5)t). (55) 

Note that rf^{t,y,x) is the fundamental solution for the operator d/dt — L*^, where L*^ is 
the formal adjoint to L^. Thus rp is the transition density for the Markov process with 
the generator Lp. Also note that ^^V'^ = 0, and thus ip'^ is the invariant density for the 
Markov process. 

Proof of Theorem \8.3[ We shall only consider the measure P/3,r since the arguments 
for the measure P f3^T{-\x{T) = 0) are completely analogous. First, let us prove the con- 
vergence of the finite-dimensional distributions. For y E M.'^ and a Borel set A G B(M.'^), 
let 

R{t,y,A) = / r[^{t,y,x)dx, 

J A 
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with rp given by (155!) . Note that R is a Markov transition function since 

r(s(t,y,x)dx = 1. 



The generator of the corresponding Markov process is Ljs and the invariant density is ip^. 
Let < ti < ... < tn < s. The density of the random vector {y'^{ti), ...,y'^{tn)) with 
respect to the Lebesgue measure on M'^" is equal to 

P^(Xi, ...,Xn) = 

Pf3{S(T) + tl, 0, Xi)pp{t2 - ti,Xi,X2)...p/3{tn - Xn-1, X„)Z;3,T-i„ (a^n) (^/3,t(0))"^ 

We replace here all factors p/j, except the first one, by r/3 using (!55l) . We replace the first 
factor and the factors Z by their asymptotic expansions given in Theorems 18.11 and 17.21 
respectively. This leads to 



p^{xi, Xn) = ^p{xi)rf3{t2 - ti, xi, a;2)...r^(t„ - + o(l), T 



00, 



where the remainder tends to zero uniformly in (xi, By the remark made after 
the statement of the theorem, this justifies the convergence of the finite dimensional 
distributions of y'^ to those of the Markov process. It remains to justify the tightness of 
the family of measures induced by the processes y'^. 

From the convergence of the one-dimensional distributions it follows that for any rj > 
there is a > such that 

P/3,t(|/(0)| >a) <r/. (56) 
for all sufficiently large T. For a continuous function x : [0,T] M"', x(0) = 0, let 

m'^(x,5) = sup |x(ti) — x(t2)|- 

\ti~t2\<5, S{T)<h,t2<S{T)+s 

Let us prove that for each e, 77 > 0, there is 5 > such that 

Pf3,T{m^{x,5) > e) <ri (57) 
for all sufficiently large T. Observe that 

Ff3,T{m^{x,S) > e) = 

rS{T)+s 

(^/3,T(0))~^Eo,T(exp( / f3v{x{t))dt)X{mT(x,5)>e}Zf3,T-S{Tysix{S{T) + s))) < 

Jo 

rS{T)+s 

(^/3,t(0))"^ sup Z^,r-5(r)-s(a;)Eo,r(exp( / f3v{x{t))dt)x{mT(x,5)>e}) < 

xSM"* Jo 

.S{T) 

exp(s/3 sup v{x)){Zp^Tifi))~^ sup Z^,T-5(T)-s(a;)Eo,r(exp( / pv{x{t))dt)x{mT{x,5)>e}) 

zeM'' Jo 
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< exp(s/5 sup v{x)){Z(s^t{0)) ^ sup Zp^T^siT)-s{x) sup pp{S{T),0,x)C{S,e), 

xeK'* xeR'i xeR'i 

where C{6,e) is the probabihty that for a c?- dimensional Brownian motion Wt, <t < s, 
we have 

sup \Wih)-Wit2)\> 6. 

\tl-t2\<S, 0<ti,t2<S 

Note that 

exp(s/3 sup v{x)){Z,3^T{0)y^ sup Zis^t-s{t)~s{x) snp Pf3{S{T),0,x) 

xeM.<i xeR'i x£R'i 

is bounded, as follows from Theorems 17.21 and 18.11 while C{6,e) can be made arbitrarily 
small by selecting a sufficiently small 6. This justifies (!57|) . Since the inequalities (156|) and 
(!57|) hold for all sufficiently large T, by choosing different a and 5, we can make sure that 
they hold for all T. Thus the family of measures induced by the processes is tight. □ 

Remark. If instead of (15^ we assume that S{T) = 0, the result of Theorem 18.31 will hold 
with the only difference that the initial distribution for the limiting Markov process will 
now be concentrated at zero, instead of being the invariant distribution. 

9 Behavior of the Polymer for p < per 

First, we shall study the asymptotic behavior of the solution Uj3{t, x) of the Cauchy prob- 
lem and of the fundamental solution pfs{t,y,x) when t — > cxd, \y\ < e~^, e\/t < \x\ < 
e~^\/i, and e > is small but fixed. Recall that (p/s was defined before Theorem 17.21 

Lemma 9.1. Let d > 3, < (3 < per, e > and f G CexplM'^), / > 0. We have the 
following asymptotics for the solution up of the Cauchy problem with the initial data f : 

up{t,x) = {2Tit)-'''^eM-\x\^/2t){{l + ^pJ)L2(^^,) + qf{t,x)), (58) 

where for some constant Cp{e) we have 

sup \qf{t,x)\ < C;3(£)r^/'||/||a.^(M3), t > 1. 
eVt<\x\<e-'^Vt 

We have the following asymptotics for the fundamental solution of the parabolic equa- 
tion: 

ppit,y,x) = (27rt)-^/2gxp(-|x|V2t)(l + yp^(y) + g(t,|/,x)), (59) 

where 

lim sup \q(t, y, x\ = 0. 

*^°°|S/|<E-S ev^<|a;|<£-iv^ 
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Proof. Note that (159|) follows from (158!) since the fundamental solution at time t is equal 
to the solution with the initial data p/sit, y, 6) evaluated at time t — S (the same argument 
was used in the proof of Theorem IS.ip . Therefore it is sufficient to prove (|58l) . 

For the sake of transparency of exposition, we shall consider only the case d = 3. From 
Lemma 15.81 it follows that we can put a = in (H3l) when /3 < Per- Thus using ( l2Til and 
the explicit formula for Ro{X), we obtain 



-1 r . Iff. p~-V2\\x-y\ 

u^{t,x) = — e'\Rp{\)f){x)dX = — / e''— -g{X,y)dydX, (60) 



where 

g{\) = iI + PA{\)r'f. (61) 

By Lemma 15.11 A{X) is an entire function of ^/X. By the Analytic Fredholm Theorem, 
(/ + (3A{X))^^ is a meromorphic function of -\/A, since ^(A) tends to zero as A — +oo, 
Im(A) = 0. It does not have a pole at zero as follows from Lemma 15.41 and Remark 1 
following Lemma 15.51 Therefore, by the Taylor formula, for all sufficiently small |A|, 
A G r(0), and some c > 0, we have 

9{X) = 90 + gi{X), ||^7i(A)||a, < cv^||/||c...,(M3), (62) 



where Qq = {I+(3A{0))-^ f . Since ||(/+/3v4(A))^i||c,,p(K3) is bounded on r(0), formula ^ 
is valid for all A G r(0), but not only in a neighborhood of zero. 
Let u^p\x) be given by fl60|) with g replaced by gi. Then 

, . Iff ^-V2\\x-y\ 

uyit,x) = — / e^*— -gi{X,y)dydX+ 

P 2m ir(o) J\y\<eVi/2 27r|a; - y\ 



Iff ^~V2X\x-y\ 

TT- / / e^'— -giiX, y)dydX = h + h- 

27^2 Jno) J\y\>eVt/2 27r|x - y\ 



We change the variable Xt = ( and use the estimate l/\x — y\ < 2/{e\/t) in Ji. This 
imphes 



'r(O) J\y\<eVt/2 



27r%t2 

C^(^)II/IIWM3) 



-\/t < \x\ < £~^Vi. 



In I2 we change the variables Xt = (, x = y/iz, y = \ftu and use the estimate < 
g-(e*/2) _ This leads to the exponential decay of I/2I as t 00. Hence 

Iff ^-V2\\x-y\ 

Uf3{t,x) = — / e-^*— rgo{y)dXdy + ri{t,x), (63) 

27r«JK3Jr(o) 27T\x-y\ 
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where the remainder ri (t, x) satisfies 

sup = ||/lbexp(K3)0(t~^) as t^oo. (64) 

£Vi<\x\<£-'^^ 

The integral over r(0) in (1631) can be evaluated, and we obtain 

1 f 

Mt^^) = (27rt)3/2 J^/ go{y)dy + n{t,x). 
Since \\go\\a^p(R'i) < C||/||cexp(K3) for some constant C, we have 



"^pit^x) = j—^e 2^ J ^go{y)dy + r2{t,x), 

where r2 satisfies (1641) with ri replaced by r2. In order to prove (I58p . it remains to show 
that 

go{x)dx = (1 + (pi3{x))f{x)dx. (65) 

Since (/ + /5fi?o(0))fi'o = we have go = f — PvRo{0)go. Recall that 9?^ is the solution 
of (US]) with / = -pv. Thus 

go{x)dx= I f{x)dx+ [ [l-A<^[^ + l3vipf3]Ro{0)godx. 

Since y^/j, i?o(0)5'o = 0(l/|x|) and their derivatives are of order 0(|x|^^) as |x| 00, the 
Green formula implies 

/ l-AippRo{0)godx = [ yDf3^ARQ{0)godx = [ ippgodx. 

JR3 Z J^3 Z J^3 



Hence 



go{x)dx = / f{x)dx+ / (fpil + (3vRoiO))godx, 



which implies ( 1651 ) □ 
Next, let us study the distribution of the polymer with respect to the measure P^^t 
as T ^ 00. Consider the process y^{t) = x{tT)/\/T, <t <1. 

Theorem 9.2. Let d > 3 and < (3 < (3cr- With respect to Pp^x, the distribution of the 
process y^{t) converges as T 00, weakly in the space C([0, 1],]R'^), to the distribution of 
the d- dimensional Brownian motion. With respect to V p;t{-\x{T) = 0), the distribution of 
the process y'^{t) converges as T 00, weakly in the space C([0, 1], M*^), to the distribution 
of the d-dimensional Brownian bridge. 
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Proof. We shall only prove the first statement since the proof of the second one is com- 
pletely similar. First, let us prove the convergence of the finite-dimensional distributions. 
Clearly P/3^t(i/"^(0) = 0) = 1. Let < ti < ... < t„ < 1. The density of the random vector 
(y'^iti), y^{tn)) with respect to the Lebesgue measure on M'^" is equal to 

T'^PpihT, 0, xiTi)p^((t2-ti)T, xiT5, X2T5)...p^((t„-t„_i)T, a;„_iT^ x„Ti)(Z^,T(0))-^ 
By Lemma [9.11 

Pp{hT, 0, xiT^) = T-'"\2Tih)-'"\l + (^^(0)) exp(-|xi| V2ti)(l + r(T, xi)), 

where 

lim sup (|r(r, xil) = 0. (66) 

Note that > po since v is non-negative, and limT^oo(-^/3,T(0)) = (1 + (/'/^(O)) by Theo- 
rem [7]2l Therefore, 

p^(Xi, ...,Xn) > 

(27rti)~5e- — (l + r(T,xi))(27r(t2-ti))~2e~Wi^...(27r(t„-t„_i))"5e 2{t„-t„_i) (57) 

= P^,...,tA^u ■■■,Xn)il + r{T,xi)), 

where t„i^iy ■■■,Xn) is the density of the Gaussian vector (W(ti), W{tn)), where W 
is a (i-dimensional Brownian motion, and q{T,xi) satisfies ( !66l) with g instead of r. Since 
e was an arbitrary positive number, this implies the convergence of the finite-dimensional 
distributions of to the finite-dimensional distributions of the Brownian motion. Indeed, 
the estimate from below for p^{xi, x„) in (167|) is sufficient since we know a priori that 
PtT... t„(^i) ^n) is the density of a probability measure. 

It remains to prove tightness of the family of processes y-^, T > 1. 

For a continuous function x : [0, T] — M'^, let 

m(x,5)= sup — x(t2)|/V^, 

|ti-t2|<<5T, 0<ii,t2<T 

m(a;, 5, £)= sup \x{ti) — x{t2)\ / \/T . 

\ti-t2\<ST, 0<ti,t2<T, |x(ti)|>ev^ 

The tightness will follow if we show that for each e,ri > there is 5 > such that 

Pi3,T{m{x,6) > e) <r] 

for all sufficiently large T. Note that m(x, 5) > e implies that m(a;, 5,5/4) > e/A. There- 
fore, it is sufficient to show that 

Vp^T{m{x,5,e/A)>s/A)<'n. (68) 
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Fix e > 0. For a continuous function x : [0, T] W^, let 

T = min(T,inf{t > : \x{t)\ = eVT/A}), 

Let £s be the event that m{x,S) > e/4 and the event that m(x, 5, e/4) > e/4. For 
< s < T, let £| be the event that a continuous function x : [0, T — s] — >• M'^ satisfies 



sup \x{ti) - x{t2)\/VT > e/A. 

\ti-t2\<ST, 0<ti,t2<T-s 



Then 



PfsA^s) = (Z/3,T(0))-^Eo,T(exp( / I3vixit))dt)xg^ 

Jo 

(Z^,t(0))-^Eo,t (^exp(^ Pv{xit))dt)E^J;p_^{x£j exp(^ , 



< 



Pv{x{t))dt)) 

where Eg denotes the expectation with respect to the measure induced by the Brownian 
motion starting at the point x. Since 



Eo,Texp(/ f3v{xit))dt) < Zfs^riO) 
Jo 

and 

Eo?L(x^jexp(/ I3v{x{t))dt)) < sup Eg,r(x^, exp( / ^v{x{t))dt)), 

Jo xm.'^,\x\=eVT/A Jo 

it is sufficient to estimate 

sup Elj.{x£,exp{[ I3v{x{t))dt)). (69) 

a;eR'*,|a:|=£v^/4 -'0 

Let 8' be the event that a trajectory starting at x reaches the support of v before time 
T. Note that 

lim sup PS,t(^') = 
since c? > 3. The expression in (1^9]) is estimated form above by 

sup {ElriXe' exp( / (3v{x{t))dt)) + PS,t(^5))- 

The second term does not depend on T due to the scaling invariance of the Brownian 
motion, and can be made arbitrarily small by selecting a sufficiently small 5. Due to the 
Markov property of the Brownian motion, the first term is estimated from above by 

sup PS,t(^') ■ sup Zp^t{x), 

xeM'',|x|=£T/4 ' xesupp(D) 

and thus tends to zero when T — > cx). □ 
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10 Behavior of the Polymer for p = per 

In this section we assume that d = 3. Again, we start with the asymptotic behavior of 
the solution up{t,x) of the Cauchy problem and of the fundamental solution pi3(t,y,x) 
when t — > oo, \y\ < e~'^, e^/t < \x\ < e^^^/t, and e > is small but fixed. 

Recall that ip is the positive ground state for Hp^^. normalized by the condition 
I l/^cr'^'i/'l (r3) = 1 (see the remark following Lemma 15.41 and Theorem 17. 2p . For / G 
C,^p(M3), define 

"(/) = ^ / 'ip{x)f{x)dx, x = 



\/2n(3cr /]g3 v{x)'ip{x)dx 

We can formally apply this to / being the (5-function centered at a point y, and thus 
define 

a{6y{x)) = xi){y). 

Theorem IQ.l. Let d = ?,, (3 = (3„, e > ^ and f G Cg^plM^), / > 0. We have the 
following asymptotics for the solution up of the Cauchy problem with the initial data f : 

upit,x) = ^-=exp(-|x|V2t)(a(/) + g;(t,x)), (70) 
\x\Vt 

where for some constant Cp{e) we have 

sup \qf{t,x)\ < C;3(£)r'/'||/||a.p(R3), t > 1. 

We have the following asymptotics for the fundamental solution of the parabolic equa- 
tion: 

Pl3{t,y,x) = —^exp{-\x\'^/2t){^{y) + q{t,y,x)), (71) 

\x\Vt 



where 



lim sup \q{t, y, x\ = 0. 

*^°°|S/|<£-S eVi<\x\<s~'^Vi 



Proof. As in Lemma 19.11 formula flTT]) follows from (1701) . Lemma 15.71 implies 

(/ + P.rA{X)r' = B + 0(1), A 0, A G C, 

P[a) — Per 

where B is the one dimensional operator with the kernel 

v{x)ij{x)ilj{y) 
}^3v[x)i)\x)dx 

From here, (l33ll and (138!) we get 

(/ + /5e.A(A))-^ = —^B + 0(1), A ^ 0, A G C, 

Pcr7vA 
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where 7 is defined in (I39]), (gO]). Hence, for any / G Ccxp(M3) and A ^ 0, A G C, 

MA, X) := (/ + ^e.A(A))- V = ^v{xMx) + ,,(A), 5(/) = 

VA /3cr(jK3^^(a;)V^(a;)dx)^ 

(72) 

where 5'i(A) < c| |/| |(7^^p(K3) for some constant c. Now, similarly to fl60|) . we have 



_i 1 /■ /■ , p-V2X\x-y\ 

up{t,x) = — / e'\Rp{X)f){x)dX = — / e^*— M\y)dydX. 

2m yp(o) 27rz Jr(o) Ju^ ^^F - Z/l 

The integral with gi{X) instead of h can be estimated similarly to the estimate on u^^^ in 
the case of /3 < Per- This leads to following analogue of (1^ 



a(f) f f Q~V2X\x-y\ 

up{t,x) = ^ / e^*— -v{y)^^{y)dXdy + n{t,x), 

2t^1' Jr3 Jy{o) 2'k\/X\x - y\ 

where the remainder ri (t, x) satisfies 



sup I 

ev^<|x|<£-iv^ 



ri{t,x)\ = ||/||a.,(M3)0(r3/2) as t 00. (73) 



We evaluate the integral over r(0): 

/ = dX=^=e ~. (74) 

2TTi ir(o) V2A 

This equality simply means that the inverse Laplace transform of the Green function of 
the one dimensional Helmholtz equation coincides with the fundamental solution of the 
corresponding heat equation. Thus, 

a(f) f 1 

Ui3{t,x) = / -e —v{y)i:{y)dy + ri{t,x). 

27r3/Vt Jr^ \x-y\ 



This implies f jTOj) since v has a compact support. □ 
The next theorem concerns the fundamental solution when both y and X are at a 
distance of order \/i away from the origin. Note that now there are two terms in the 
asymptotic expansion for the fundamental solution which are of the same order in t. The 
main terms have the order when t ^ 00, compared with in the case considered 
in Theorem 110.11 (where y was bounded). 

Theorem 10.2. Let d = 3, (3 = jScr, £ > 0. We have the following asymptotics for the 
fundamental solution of the parabolic equation: 

Pfsit, y, X) = po{t, y, X) + } ,-MMA?/2t^^ + ^^^^ ^))^ (75) 

[2TiYi^\y\\x\yt 



where 

lim 

eVt<\y\,\x\<e-^yft 



lim sup |g(t,y,x| = 0. (76) 

t^oo 
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Proof. Let Pjsit, y, x) = po{t, y, x) + u. Then Ut = Hpu + /?fpo, ^14=0 = 0, and therefore 
by the Duhamel formula 

u{t,y,x)= / pf^it - s,z,x)(3v{z)po{s,y,z)dzds. 



Usine; fITT]). we get 







u{t,y,x) = I I =exp{—\x\'^ /2{t — s)){ilj{z)j3v{z)po{s,y,0)dzds + hi + h2 (77) 



Jr^ - s 



with 



hi= / I I ^ exp{-\x\'^/2{t - s)){tlj{z)(3v{z){po{s,y,z) - po{s,y,0))dzds, 

Jo \x\y/t - S 

h2 = / =exp{—\x\'^/2{t — s)){q{t — s,z,x)f3v{z)pQ{s,y,z)dzds, 

Jo Jr3 mvt - s 

where q is the same as in (!7T|) . The integral in the right hand side of (177|) (let us denote it 
by w) is a convolution of two functions and can be evaluated using the Laplace transform 
(see ( FM|) ). It gives the second term in the right hand side of (175|) . The contribution from 
the other two terms can be shown to satisfy (1761) . Let us prove the statement about w. 
In fact, 



XV 27r f 11 

w = >ci{wi*pQ(t,y,0)), xi = — j — j — / l3v{z)ip{z)dz = Wi = —^=exp(—\x\'^ /2t)). 

fI Jr3 \x\ V27rt 



The Laplace transform Wi{X) of the function wi is equal to e ^^l^'l/v^2A (see fITil) ). and 
the Laplace transform of pQ{t,y,0) is equal to e~^^'^'/27r|?/|. Thus 



1 V2A{|x|+|s/|) 

^;(A) 



It remains to apply (17^ one more time. □ 

As in Section [9l we shall study the limit, as T —>■ oo, of the family of processes 
yT(^t) = x{tT)/VT, 0<t<l. For < s < t < 1, ?/,a; G M^ define 

pj(s, t, y, x) = Pp{T{t - s),yVT, xVT), 

Pfsis, t, 0, x) = lim (TpJ(s, t, 0, x)) = lim {Tpp{T{t - s), 0, xv^T)), x ^ 0, 

pJs, t, y, x) = lim {T^''^pl{s, t, y,x)) = lim iT'^/'^pp{T{t - s),yVf, xVt)), y,x^O. 
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By Theorems [TO and [TO 

Pf,{s,t,0,x) = ^^^2Lexp(-|x|V2(t - s)), x^O, 
\x\\'t — s 

Pp{s, t, y, x) = po{t - s, y, x) + ] exp(-(|y| + \x\)^/2{t - s)) ,y,Xy^ 0. 

[27Tj'^'^\y\\x\\'t — s 

(78) 

For < ti < ... < tn < 1, let the density of the random vector {y'^{ti), ...,y'^{tn)) with 
respect to the Lebesgue measure on M*^" be denoted by p^(xi, ...,x„). 
For < s < t < 1 and y,x E M^, define 



Q^{s,t,y,x) =p'^{s,t,y,x) p'J^{t,l,x, z)dz{ pl{s,l,y,z)dz) ^, t < 1, 

0-(M....)^P?(M....)(/pJ(M....M.)- 



Thus 

P^(Xi, ...,Xn) = (5^(0, ti, 0, Xi)Q'^(ti, t2, Xi, X2)--Q^{tn~l,tn, Xn-l, Xn). 

In order to find the hmit of the finite dimensional distributions of y'^ , we need to identify 
the limit of Q"^ as T ^ oo. For < s < t < 1, ?/ G and x G \ {0}, define 

Q(s, t, y, x) = lim Q^(s, t, y, x). (79) 

T—>oo 

By Theorems [TO and [TO 

(81) 

We additionally define Q{s,t,y,0) = 0. 

Using (17^ . (I5U1) and (IHTl) . we can identify the limit of the densities p^{xi, ...,x„) for 
X2, ■■■,Xn 7^ 0. In order to identify the weak limit of the finite dimensional distributions 
of the processes y'^ , we are going to show that the limit of the densities is the density of 
a probability distribution, i.e. the mass does not escape to the origin or infinity. This 
is done in Lemma 110.41 where we show that Q serves as the transition density for a 
Markov process. First, however, we show that Q satisfies a Fokker-Plank type equation 
on \ {0}. 

Let 

g{t,x) =\ii{ pp{t,l,x,z)dz), 0<t<l, |x| > 0, (82) 

JRS 
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Let L be the differential operator acting on C^(]R^ \ {0}) according to the formula 

(L/)(t,x) = \^J{t,x) + |x| > 0, 

and let L* be the formal adjoint of L, i.e. 

1a 1 d\{dg/dr)v\ 



2 dr 

Lemma 10.3. For < s < 1 and y G M? , the function Q{s,t,y,x) satisfies the equation 

^9i^2liyi^ = L*Q{s,t,y,x), \x\>0, s<t<l. (83) 
Proof. Let us consider the case when y ^ (the other case is similar). Let 

vi{s,t,y,x) = I exp(-(|y| + \x\f /2{t - s)), 

{27Ty/^\y\\x\Vt- s 

V2it,x)= [ ] ^ =expi-{\x\ + \z\y/2{l-t))dz. 

Jr3 (27r)'^/^|a;||z|Vl - 1 

Observe that 

(|-1a>. = 0. (| + iAJ„. = 0. (84) 
For fixed s and y, the function Q{s,t,y,x) is proportional to 

u{t,x) = {po{t- s,y,x) +Vi{s,t,y,x))[l + V2it,x)]. 

By (El, 

For any two functions A and B we have 

(a|- + B)u = + + ^2) + A{po + v,)^ + Bipo + v,){l + V2). (86) 

or or or or 

Thus 

5 1 d 



if 



Since g{t,x) = ln(l + ^2) and 2dv2/dt = — 5^f2/5r^ — 2dv2/dr (see fl84l)). it is easy to 
check that the operator in the left hand side of the equation for u is — L*, and this 
justifies dHS]). □ 
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Lemma 10.4. The function Q{s,t,y,x), < s < t < 1, y,x & M^, is the transition 
density for a Markov process on M.^ . 

Proof. To show the existence of a Markov process, we need to verify that 

/ Qiti,t2,Xi,X2)dX2 = I, ti<t2 (87) 

and 

Q(ti,t2, Xi, X2)Qit2, ts, X2, X3)dx2 = Q{ti, ts, Xi, X3), ti < ^2 < ^3- (88) 



Let us assume that ( !87l) has been demonstrated, and prove ( l88l) . Observe that 

/ T^+>^(ti, t2, Xi,X2)p]s{t2, ts, X2, X3)dX2 = T^+>^()f:i, ta, Xi, X3), ti < t2 < t3, 

where a = 1/2 if xi = and a = otherwise. For 0:3 7^ we take the hmit, as T ^ 00, 
on both sides of this relation. The integrand on the left hand side converges to 

Pp{ti,t2, Xi, X2)pp{t2, ts, X2, X3), 

however, the convergence is not necessarily uniform in X2, and we can only conclude by 
the Fatou Lemma that 

Pp{ti,t2, Xi, X2)Pfi{t2, ts, X2, X3)dx2 < P/3(tl, ^3, ^^i, X3), ti < < ^3, X3 7^ 0. 

IR3\{0} 

From ([80]) and ^ it now follows that 

Q{tl,t2,Xi,X2)Q{t2,t2.,X2,X3)dX2<Q{ti,t3,Xi,X3), ti < t2 < ^3, X3 ^ ^. 



Note that both sides of this inequality are continuous in X3 e \ {0}. Due to (IHTI) . the 
integrals in X3 over \ {0} are equal to one for the expressions in both sides of this 
inequality. Therefore, 

/ Q{ti, t2, Xi, X2)Q{t2, t3, X2, X3)dX2 = Q{tl, ^3, ^^l, X3), ti < ^2 < ^3, X3 ^ 0, 

Jr^\{o} 

and thus (EZD implies (l88l). 

Now let us verify ( IHTl) . Put s = ti,T = t2,y = xi and x = X2. Again, we shall consider 
the case y ^ 0, the other case being similar. Moreover, we can assume that t < 1, since 
the case r = 1 can be treated by taking the limit t | 1. On a formal level, (IHTI) follows 



from (l83l) by integrating the both sides of (183|) over f2 = [s, r] x C 

/ Q{s, T, y, x)dx -lim Q{s,t,y,x)dx = {L*Q,1)l2^<^^ = {Q,L1)l2(^<^). (89) 
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One needs only to note that 



lim / Q{s,t,y, x)dx = 1, (90) 

and that the operator L apphed to the identity function gives zero. The latter implies 
that the left hand side in fl89|) is zero, and (l90l) implies that the second term on the left 
hand side of fl89l) is one. 

In order to make relations fl5^ rigorous we note that Q{s,t,y,x) is infinitely smooth 
in {t,x) when x and decays exponentially as |a;| — > oo. However, it has a singularity 
at a; = 0. Thus the integrals over and Q in (jHlD must be understood as limits of the 
corresponding integrals over the region |x| > e as e — >■ 0. Let us examine the singularities 
of Q and of the coefficients of L* at the origin. 

Relation (ffSll implies that 

Pf^{s,t,y,x) = ^ + 0{r), r = |x| ^ 0, a = a{s,t,y). (91) 

It is important that (l9T|) does not contain a term of order 0(1). From (I9TI) . (!82l) and (l80l) 
it follows that 



dg(t,x) 1 ^ ^/ i „\ ^ , 9Q{s,t,y,x) 2c 

dr r 



3 + 0(r), Q(.,t,y,x) = ^ + 0(1), ^ =-^ + 0(1), (92) 



where r — 0, c = c{s,t,y). Since Q has a weak singularity at x = 0, the integral of the 
left hand side of (1831) over = Q f]{x : |x| > e} converges to the left hand side of (l89l) . 
Hence, in order to prove (|87ll . it remains to show that 

L*Qdtdx ^0, e^O. 

The integral above is equal to 

where da is the element of the surface area of the sphere |x| = e. The convergence of fl93l) 
to zero follows immediately from ( 192|) □ 

Lemma 10.5. The family of processes y'^it), T >\, is tight. 

We shall prove this lemma below. First, however, we formulate the main result of this 
section. 

Theorem 10.6. The distributions of the processes y'^{t) converge as T oo, weakly 
in the space C([0, 1], M^), to the distribution of the ?>- dimensional Markov process with 
continuous trajectories. The transition densities for the limiting Markov process are given 
by (Eg) and 
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Proof. The convergence of the finite dimensional distributions of y'^{t) to those of the 
Markov process follows from (!79l) and Lemma [10. 4[ Since the family y^{t) is tight, there 
is a modification of the Markov process which has continuous trajectories. □ 

Proof of Lemma \10.5[ To prove tightness it is enough to demonstrate that for each 
ri,e > there are < 5 < 1 and Tq > 1 such that for all u G [0, 1] we have 

F/,,t{ sup \y^{s)-y^{u)\>e)<6r], T > Tq. (94) 

ti<s<min(t+5,l) 

Let ri,6 > he fixed. Let Ss be the event that a continuous function x : [0, T] M.^ 
satisfies 

sup \x{t) - x{0)\/Vt > e/8. 

t<ST, 

Using arguments similar to those leading to (l69l) . we can show that (p4|) follows from 

sup Ey{x£,expi[ Pvix{t))dt)) < Sri, T > Tq. (95) 

x£R'i,\x\=eVT/4: JO 

Let 

r = min(5r,inf{t > : |x(t) - x(0)| = eVT/S}), 
The expectation in ( 195|) can be estimated as follows 

pT pT—t 

KAxes exp( / f3vi4t))dt)) < ES,r(x^,E;^5L. exp( / /3t;(x(t))rft)) 
Jo Jo 

We claim that 

EofL^exp(/ (3v{x{t))dt) < sup El^expi (3v{x{t))dt) < c{e) (96) 

"'0 xm'i,\x\>eVT/8 ' Jo 

for some constant c{6) for all sufficiently large T. It then remains to choose 6 such that 
^orixsg) < '^^/c(^), and the estimate (195!) will follow. The second inequality in (!96l) 
easily follows from part (2) of Theorem 17.21 and the fact that the probability of reaching 
the support of v before time T by a Brownian path starting at a distance eVt/8 away 
from the origin is of order 0(T~^/^) ii d = 3. □ 
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